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We introduce a set of discrete modular transformations Ti, Ui and St in order to study the relation- 
ships between the different phases of the Heisenberg ladders obtained with all possible exchange 
coupling constants. For the 2 legged ladder we show that the RVB phase is invariant under the 
St transformation, while the Haldane phase is invariant under Ut- These two phases are related by 
Ti. Moreover there is a "mixed" phase, that is invariant under T(, and which under Ut becomes 
the RVB phase, while under St becomes the Haldane phase. For odd ladders there exists only the 
Tt transformation which, for strong coupling, maps the effective antiferromagnetic spin 1/2 chain 
into the spin 3/2 chain. 

PACS numbers: 75.10.Jm, 75.50.Ee 
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In the last two years the concept of duality has played 
a crucial role in understanding non perturbative_aspects 
of Quantum Field Theorytl and String Theorytl. Con- 
sidering the traditional links between Particle Physics 
and Statistical Mechanics or Condensed Matter one may 
wonder wether these latter areas could benefit from the 
deeper understanding gained in the former ones. In 
fact, duality ideas have been important in the histori- 
cal development of Statistical Mechanics, as shown by 
the Krammers-Wannier duality, order-disorder transfor- 
mations, eto9. In this letter we shall explore the existence 
of duality symmetries in quantum spin systems defined 
on a lattice and more particularly on arrays of coupled 
spin chains known as spin ladders! 

Generally speaking a duality transformation is a map- 
ping between two models, or the same model with differ- 
ent parameters, which apparently have different physical 
properties, but which become in a way equivalent under 
the transformation. Dual theories usually give comple- 
mentary descriptions of the same underlying phenomena. 

Let us first stablish what we mean by duality in a spin 
system. We shall consider the Heisenberg Hamiltonian 
defined on the <i-dimensional hypercubic lattice (d > 1), 



#(W) = E E j a< s *- s > 



(i) 



where S x is a spin S matrix acting at the position x = 
(xi, . . .,Xd), and^i = (1, 0, . . . , 0), . . . , n d = (0, 0, . . . , 1). 

We shall define the dual of the Hamiltonian (Q) as a 
Hamiltonian Ho = H({J^}) characterized by a new set 
of coupling constants { }, and such that the low energy 
spectrum of H and Ho is in one to one correspondence. 
This implies that the free energy and the ground state 
energy will also be the same for both models. 

In the classical limit where the spin S> 1, the ground 
state of (Q) is given by the classical minima, 



(2) 



where n is a 3-component unit vector and e M = 
The energy of ^ is given by, 



-sign J... 



= -S* 
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(3) 



The signs of the exchange coupling constants, e M , de- 
termine the type of order parameter which characterizes 
the ground state. Thus e p = 1 or —1 correspond to 
ferromagnetic (F) or antiferromagnetic (A) order in the 
direction /j, of the lattice. Altogether there are 2 d pos- 
sible classical vacua which we denote, for d — 2, by the 
sequence, 



AA, AF. FA, FF 



(4) 



The energies of the classical g.s. and the excitations 
are independent of the type of vacua (^). All the clas- 
sical Heisenberg models are equivalent. However the 
quantum corrections drive them into very different quan- 
tum vacua. Only the pure ferromagnetic system (i.e. 
J M < 0,V/i), survives the quantum fluctuations, but the 
non- ferromagnetic systems change deeply their structure. 
The purpose of this letter is to show the relations existing 
between the different vacua by means of a certain type 
of duality transformations. 

At this point it is useful to make an analogyJpetween 2d 
spin systems and fcrmions living on a 2d torusO. To define 
a fermion on a torus one has to specify the boundary 
conditions along the a and b cycles. They can be periodic 
(P) or antiperiodic (A) . This gives rise to 4 possible spin 
structures, labeled as AA, AP, PA and PP, which mix 
under the actios, of the modular transformations T, U 
and S as followsa, 



T : AA «-> AP, PA, PP 
U : AA^ PA, AP, PP 
S : AP ^ PA, AA. PP 



(5) 



Observe that the spin structure PP is left invariant un- 
der the action of the modular group. The fermion deter- 
minant with the boundary conditions AA, AP, PA turns 
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out to be given by Jacobi functions, which transform 
among themselves under the modular group as described 
by (|]). The fermion determinant for PP boundary con- 
ditions is zero due to the existence of a zero mode. 

In the case of 2d spin systems the role of the cycles a 
and b is played by the directions ji\ — (1,0) and /i 2 = 
(0, 1). The analog of the spin structure is given by the 
(anti) ferromagnetic nature of the coupling J M along the 
directions ^1,2- Finally, a modular transformation is a 
redefinition of the unit cell of the lattice. In the case of 
spin ladders the above analogies can be collected in the 
following dictionary, 



Torus Lattice 
a — cycle 
b — cycle 
Antipcriodic BC 
Periodic BC 
Modular Transf . 



Spin Ladder 

legs 

rungs 

Antifcrromagnetic Coup. 
Ferromagnetic Coup. 
Bond Moving Transf. 



(6) 



These correspondences have an analogue in d > 2. On 
a 2 legged ladder, 2-ladder for short, we shall define three 
transformations Ti, Ug and Si as follows. The Tg trans- 
formation consists in the shift by one lattice spacing of 
one leg respect to the other one (see fig.l). The Ug trans- 
formation consist in the permutation of the two sites of 
the even rungs, while leaving invariant the odd ones (see 
fig. 2). Finally the Si transformation is defined by the 
equation Si — TgUgTg, and has the effect of converting 
all the vertical bonds (rungs) into horizontal ones (legs), 
while half of the horizontal bonds become vertical bonds 
and the other half become diagonal bonds of length v5- 
We remark that Tt, Ug and St do not generate the stan- 
dard modular group. 

Using these definitions one can see that the classical 
vacua AA, AF, FA and FF, get mixed under the action 
of Tg,Ug,Sg in the form described by (j^), with the re- 
placement: (anti)periodic «-> (anti)ferromagnetic. The 
term "bond moving" in @ refers to a transformation 
introduced by Migdal and Kadanoff in the study of the 
Ising model with RG methods!! Just like a fermion on a 
lattice with non periodic BC's can be regarded as essen- 
tially the same object, we conjecture that the 2-ladder, 
and more generally the ladders with an even number of 
legs, with "magnetic structures" , AA, AF, FA, belong to 
the same universality class and are related through dual 
transformations. This conjecture implies in particular 
the equivalence between the RVB state (AA couplings) 
and the Haldane state (AiTpCOuplings) of the 2-ladder 
studied by various authors3,Q, but it also suggests new 
equivalences which have not been studied so far. We shall 
confine ourselves in this letter to the case of the spin 1/2 
ladder with 2 legs, trying to prove the above conjecture 
using perturbative and field theoretical techniques. At 
the end we shall briefly consider the case of odd ladders. 

Let us start with a toy ladder. 

A 2 x 2 cluster 

The simplest 2-ladder has 4 spins coupled by the 
Hamiltonian. 



H = J a (Si • S 2 + S 3 ■ S 4 ) + J 6 (Sa • S 4 + S 3 ■ S 2 ) (7) 

The ground state of the non-ferromagnetic Hamiltoni- 
ans (i.e. e a ,b 9^ 1) is a singlet and therefore can be written 
as the linear combination 



m =T\a)+ \b) 
\a) = (12)(34), \b) = (14)(32) 



(8) 



where (£7) denotes the singlet valence bond state con- 
structed out from the spin 1/2 's located at the sites i 
and j. The transformations Sg,Tg and Ug of figure 1 be- 
come for the toy ladder elementary transpositions, 



Tg : 3 <r* 4, Ug : 2 <-> 3, Sg : 2 4 



(9) 



The action of Tg,Ug,Sg on the states (||) can be easily 
derived from (@) and (ph. They are given in table 1. 





Tg 


Ug 


Sg 


\a) 
\b) 


-|a> 
-|o) + \b) 


\a) - \b) 
-\b) 


\b) 
\a) 



Table 1. 

The "modular transformations" induced on the "mod- 
ular parameter" r that follow from table 1 are, 



1), Ug-.r^ -t/(t + 1), Sg:T^ 1/t 

(10) 



which are similar but not identical to the standard mod- 
ular transformations of the torus. 
The ground state energy of (^) is, 



E 
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(Ja + Jb) - \ Na + Jb - J aJb 



corresponding to a value of r given by, 



(11) 



(12) 



The values of r obtained by changing the signs and 
strengths of the exchange coupling constants J a ,b, cover 
the whole real axis as described in table 2. 




(£a,£fc) 


AF 


FF 


FA 


AA 


T 


(-oo,-2) 


(-2,-1/2) 


(-io) 


(0,oo) 



Table 2. 

We have included the case FF which corresponds to 
an excited state, since the g.s. is a spin 2 multiplct. 
The Sg transformation ( [Tof ) leaves invariant the AA and 
FF domains, while interchanges the regions AF and FA. 
^-duality is an exact symmetry of the Hamiltonian (Q) . 
Actually, r = 1 is a fixed point of St- The Tg and Ug 
transformations are approximate symmetries in the sense 
that the Hamiltonian (Q) is not mapped into a similar 
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one with a redefinition of J a ,b- However one can see that 
t = —1/2 is a fixed point of Te, while the A A region 
r > 1 is mapped under Te into the AF region r < — 2. 
Similarly r = — 2 is a fixed point of Ue, while the 
region < r < 1 is mapped under Ue into the FA region 
-1/2 < t < 0. All these shows that eqs.(||) hold with 
some caveats for the 2x2 cluster. 

Within each domain, A A, AF and FA, we shall choose 
a representative state |r) with the property of being in- 
variant under one of the dual transformations. The state 
t = 1 can be called a RVB state since it describes the res- 
onance between two vertical and horizontal bonds. The 
state r = —2 is a Haldane like state (HAL) in the sense 
that it is obtained upon forming the spin 1 state along 
the rungs, which then couple to form a singlet. Finally 
t = — 1/2 is a mixed state (MIX), corresponding to fer- 
romagnetic chains coupled antiferromagnetically. More- 
over, each of the states \RVB), \HAL) and \MIX) gets 
transformed into one another by the action of Te, Ue, Se- 

The results are summarized in table 3. 



Te ■ (g,g) -> (-9,9) 



(15) 



(Ja, Jb) 


State 


T 


Se 


Ue 


Te 


AA 


RVB 


1 


RVB 


MIX 


HAL 


AF 


HAL 


-2 


MIX 


HAL 


RVB 


FA 


MIX 


-1/2 


HAL 


RVB 


MIX 



Table 3. 

For 2-ladders with a large number of rungs we can still 
make sense of the transformation properties collected in 
table 3. In that case \RVB) denotes the ground state of 
a ladder with AA couplings, etc. 

The rest of the letter will be devoted to show the va- 
lidity of table 3. 

The Weak Coupling Regime: T^-duality 

If the two legs are weakly coupled (i.e. \J a /Jb\ » 
1), the T^-duality becomes a manifest symmetry of the 
effective low energy theory. 

For J a > we can use bosonization techniques to show 
this fact. Indeed the effective ladder Hamiltonian can be 
written in the bosonized model asE3, 



H = H 



wzw 



Ai(Jj 



•J*) 



+A 2 (Ji 



\ 3 Tr(g<T)Tr(ga) + A 4 TrgTrg 

(13) 



where g and J (resp. g and J) are the WZW field and 
current which bosonize the upper (lower) spin-chains of 
the ladder. The initial values of the different coupling 
constants appearing in ( |l3| ) are given by 



Ax < 0, A 2 = A 3 = J b , A 4 = 



(14) 



In the bosonized representation the translation of r 
chain by one site is equivalent to the discrete symmetry 
g — ► —g. Therefore the operator Te is realized in the 
WZW model by the map, 



implying that Te is equivalent to the following change of 
couplings, 



7 S 7 

,J a * iJ a 
Jb —Jb 



(16) 



Eq.(|l6|) illustrates the relations Te \RVB) = \HAL) and 
Te \HAL) — \RVB), which stablish the equivalence be- 
tween the RVB and Haldane states in the weak coupling 
limit. White has observedcl that the spins located in di- 
agonal positions of the 2-ladder tend to form effective 
spins 1, and using, what we call the iy transformation, 
he shows ±he equivalence between the two phases. This 
is done incl by introducing diagonal couplings in order to 
connect continously the \RVB) and \HAL) states. What 
we show in this letter is that this connection can be also 
thought as a discrete modular transformation by which 
the properties of both models can be put in one to one 
correspondence. 

If J a < both legs are in a ferromagnetic state with to- 
tal spin S to t = N/2. A weak antiferromagnetic coupling, 
Jb > 0, splits this degeneracy giving a state which, to first 
order in perturbation theory, is given by the singlet ap- 
pearing in the Clebsch-Gordan decomposition £ t ot x <Stot- 
Obviously, the latter state is invariant under a shift of one 
of the legs. Thus the state \MIX), in the weak coupling 
regime, is invariant under Te— duality, according to table 
3. 

The Strong Coupling Regime: [/^-duality 

In the strong coupling regime (i.e. \J a /Jb\ *C 1) the 
rung Hamiltonian yields the zero order approximation, 
while the leg Hamiltonian acts as a perturbation. 

The rungs, in a AF ladder, are mostly in the spin 1 
state which couple antiferromagnetically along the leg 
direction, yielding effectively a Haldane chain. The Ue 
transformation, which simply permutes the two spins on 
the even rungs, leaves invariant the corresponding Hal- 
dane state (i.e. U e \HAL) = \HAL)). 

Next we shall show using perturbation theory that the 
RVB and MIX states are exchanged by [/^-duality. 

For Jb » 1 the rungs are in a singlet state.,-The g.s. 
energy computed to 2 nd order in J a is given bylJ, 



E /N 



;Jb 



ZJl 
8 J b 



(17) 



The first excited states form a band of spin 1 magnons, 



\k) 



1 N 



e lkx \x) 



(18) 



where \x) denotes the state with singlets on all rungs ex- 
cept at the position x where it is a triplet. The dispepion 
relation ui(k) of (|l8| ) is given, to 2 nd order in J a , bjO, 

1 J 2 

uj(k) = Jb + Ja cosfc - _ -Ji (3 - cos2fc) (19) 
4 Jb 
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The action of Ue on the magnons (|1 
Ut \k) = \k + n) 



(20) 



Hence the spectrum of the RVB and MIX states, up 
to second order in J a , are exchanged under [/^-duality, 
as can be seen from the following identities satisfied by 
©and©, 



Eo(J a , Jb) — E (—J a , Jb) 

U)(k, J a , Jb) = Oj(k + 7T, - J a , J b ) 



(21) 



The Intermediate Coupling Regime: 5^-duality 

When \J a /Jb\ ~ 1, the effective theory can be obtained 
by map ping . thje ladder into the non linear sigma model 
(NLSM)MaH The values of the NLSM coupling con- 
stants are given in table 4, 





AA 


AF 


FA 


6 





4tt5 





9 


1 




1/2 


l 

s 


1 


\2\Ja\, 


1/2 



Table 4: Parameters of the 2-ladder with spin 5. 
From these eqs. we get the curious relation, 



2 2 2 

9aa — 9af + 9fa 



(22) 



Sf-duality corresponds to the permutation of vertical 
and horizontal bonds. Since on a 2-ladder there are twice 
as much horizontal bonds than vertical ones we expect a 
perfect balance between both couplings whenever 2| J a \ = 

\Jb\. In this case gAF — 9fa = 9aa/V%- The change of 
6 by AttS, when going from FA to AF, does not affect 
the physics of the problem and recalls what happens with 
duality transformation in field theoriesu. 

Extrapolating the NLSM-map away the intermediate 
couplings we still find an agreement with table 3. In the 
strong coupling regime both qaa and gFA go to the same 
asymptotical value, which agrees with the fact that the 
Ue operation maps one ground state into the other. On 
the other hand, the value gAF = 1/5 = 2 corresponds to 
the NSLM coupling of a spin chain with spin 1. 

In the weak coupling regime the NLSM map is not 
reliable, however we see from table 4, that in that limit 
9AA = 9af 9fa which agrees with the fact that Te 
interchages the RVB and HAL states. Of course in this 
limit we have two weakly coupled chains, which should 
be treated with bosonization techniques. 

Summarizing our results we can say that the Haldane 
and the mixed phases are Se— dual, while the RVB phase 
is self-dual under a 5^-transformation. 

Bond Moving Dualities 

What is the origin of the duality properties of ladders? 
In Conformal Field Theory or String Theory duality is 
an expression of modular invariance. Something of this 
sort exists also for spin systems. To show this we shall 



use a generalization of the Migdal-Kadanoff transforma- 
tions, which consist in the substitution of couplings be- 
tween NN sites by other NN or NNN couplingsEL This 
is achieved by adding a potential V to the Hamiltonian 
H, so that the new Hamiltonian H' = H + V has a g. 
s. energy ( and free energy) E' smaller than the g.s. en- 
ergy E of H, provided (V) = 0, where the v.e.v is taken 
respect to the g.s. of HU. From fig. 1 we observe that 
the Te— transformation corresponds to the bond moving 
potential, 



V T t =J2( J b Si(n)-S 2 (n+1)- J b Si(n)-S 2 (n)) (23) 



The Hamiltonian obtained adding (23) to the ladder 
Hamiltonian H(J a , Jb), is a new Hamiltonian H(J' a , J' b ), 
where J a ,b are given by 



J' b = J 6 (S 1 (n)-S 2 (n))/(Si(n)-S 2 (n 



(24) 



1)> 



The expectation values in ( j24| ) are computed with re- 
spect to the g.s. of H(J a , ,7fc).Eq.(p4l) implies sign(J^) = 
— e a eb, which indeed corresponds to a Ti— transformation. 

Similarly the bond moving transformation which cor- 
responds to Ue gives, 



J'b 



Jb (25) 
J' a = J„{S!(n) • Si(n + l))/(Si(n) ■ S 2 (n+ 1)> 



Finally the 5^-tranformation can be derived from its 
definition Se — TeJJeTe- For all these transforma- 
tions, duality would amount to the equality E(J a , Jb) = 
E(J' a ,J'b). The variational principle underlying the 
Migdal-Kadanoff transformation only guarantees that 
E(J a , Jb) > E(J' a , J' b ), however after the results obtained 
above using perturbative and field theoretical methods, 
we have good reasons to believe that the replacement of 
inequalities by equalities for the energies and free ener- 
gies yields a good approximation. Further studies are 
necessary to fully stablish these facts. 

Beyond the 2-ladder Most of the results shown so 
far are generalizable to the case of even-ladders with pe- 
riodic BC's along the rungs. The T^-transformation is 
given by the shift of one lattice space of the even legs 
respect to the odd legs, so that the rungs become zig- 
zag lines across the ladder. In fact, this definition ap- 
plies to all types of ladders, even and odd, with different 
BC's across the rungs. The Ue transformation consists in 
the shift by one lattice space of the even rungs, so that 
the legs become zig-zag lines along the ladder. The new 
coupling constants, obtained upon these transformations, 
are also given by the eqs.(p^) and j25|). 

For odd-ladders with open BC's along the rungs, there 
seems to be not a sensible definition of the Ue and Se 
transformations, as we did above for Te. The odd-ladders 
of type AA are in the same universality class as the spin 
1/2 antiferromagnetic Heisenberg chain, whose g.s. we 
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shall denote as jA^). On the other hand the n^-ladders 
(nf. odd) with AF couplings are, at least in the strong 
coupling regime, equivalent to spin ni/2 antiferromag- 
netic chain, whose g.s. we shall denote as \A ni / 2 )- The 
role of Ti is to exchange the A A and AF couplings, which 
implies that Tt I-A1/2) = \A ne /2)- The previous equiv- 
alence can be stablished for every regime of couplings 
using, as we did for the 2-ladder, the appropiated tech- 
nique. Thus for intermediate couplings, wheca we can 
use the mapping of the ladder into the NLSML3, we get 
the parameters &aa = tt and Oaf = ^ng, which coincide 
modulo 2ir. The odd-ladders with FA and FF couplings 
are equivalent to ferromagnetic Heisenberg chains with 
spins 1/2 and ne/2 respectively. Their g.s. is invariant 
under T&. 

In summary we have seen that the even and odd lad- 
ders have quite different duality properties which is of 
course a manifestation of the fact that they both belong 
to different universality classes. 
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Figure 2 



